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ABSTRACT 

I ■ Context. H2O is the most abundant component of astrophysical ices. In most lines of sight it is not possible to fit both the H2O 3 fim stretching, 
' the 6 fim bending and the 13 fim libration band intensities with a single pure H 2 spectrum. Recent Spitzer observations have revealed CO2 
f — . ice in high abundances and it has been suggested that CO2 mixed into H 2 ice can affect the positions, shapes and relative strengths of the 3 
f^*) ' fim and 6 fim bands. 

r- H Aims. We investigate whether the discrepancy in intensity between H2O bands in interstellar clouds and star forming regions can be explained 
. by CO2 mixed into the observed H 2 ice affecting the bands differently. 

' Methods. Laboratory infrared transmission spectroscopy is used to record spectra of ILO:C02 ice mixtures at astrophysically relevant 
(-H i temperatures and composition ratios. 

O I Results. The H2O peak profiles and band strengths are significantly different in H20:C02 ice mixtures compared to pure H2O ice. The ratio 
between the strengths of the 3 fim and 6 fim bands drops linearly with CO2 concentration such that it is 50% lower in a 1:1 mixture compared 
to pure H 2 ice. In all H 2 0:C02 mixtures, a strong free-OH stretching band appears around 2.73 fim, which can be used to put an upper limit 
on the CO2 concentration in the fLO ice. The fLO bending mode profile also changes drastically with CO2 concentration; the broad pure ILO 
^ band gives way to two narrow bands as the CO2 concentration is increased. This makes it crucial to constrain the environment of H2O ice to 
, enable correct assignments of other species contributing to the interstellar 6 fim absorption band. The amount of CO2 present in the H 2 ice of 
^ ' B5:IRS1 is estimated by simultaneously comparing the H 2 stretching and bending regions and the CO2 bending mode to laboratory spectra 
of H 2 0, CO2, H 2 0:C0 2 and HCOOH. 
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1. Introduction 

Infrared spectroscopy towards dense molecular clouds and 
young stellar objects (YSOs) often reveals prominent bands 
attributed to H 2 ice. The 3.08 fim (3250 cm" 1 ) band, at- 
tributed to the H2O stretching mode, and the 6.0 fim (1700 



in many lines of 


sight < Willner etal. 1982; 


Tanaka et al||l990} 


Murakawa et al. 


2000; 


Booaert et al. 


2000 


Keaneet al. 2001; 


Gibb et alJ2004llKnez et alJ2005h. Also the observation of the 



H7O libration mode at 13 ^m (770 cm -1 ) has been reported 
( Kea ne et alJ2 001: Gibb et al. 2004|). It is a long-standing prob- 
lem in the interstellar community that the H2O ice 6.0/3.08 
band intensity ratio in astrophysical observations is up to 2 



times higher than expected for pure H2O ice jKeane et alJ 
2001). 
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This discrepancy has previously been explained by absorp- 
tions of other volatile molec ules and organic refractory mate- 
rial absorbing around 6 fim dSchutte et al.lll996l IKeane etafl 
200l[ lGibb & Whittetll2002l) . These theories were put forward 
to explain the strong absorption of the 6.0 fim feature as well 
as why its shape does not match that of pure H2O ice very well 
in many lines of sight. IKeane et"al] d200ll) identified two ad- 
ditional components in the H2O 6.0 fim bending mode region, 
centered around 5.83 fim (1720 cm -1 ) and 6.2 fim (1600 cm -1 ), 
when subtracting a pure H2O ice spectrum that was fitted to 
the 3.08 fim sketching band. Analyzing the residuals after sub- 
tracting the H2O bending mode presupposes that the profile of 
the H2O bending mode is well known. Hence knowledge of the 
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H2O bending mode profile in different environments is critical 
to correctly assign other species contributing to the 6 fim band. 

IPontoppidan eta showed that for observations of 

ices in circumstellar disks, part of the discrepancy in intensities 
between the H2O bands may be due to disk geometry because 
of scattering at short wavelengths. However, since this anomaly 
is almost always present to some degree in disks and envelopes, 
as well as in clouds, it is unlikely to be the entire explanation. 
In the sources where the 13 fim H2O libration band is visible, 
it is possible to fit the H2O stretching and libration peaks with 
a single H2O ice abundance, while fitting the 6 fim band inde- 
pende ntly results in a much higher column density ( Gib b et all 
12004 . The fact that the libration mode at 1 3 fim is in agreement 
with the 3 fim band refutes the idea that the excess at 6 fim is 
due to wavelength dependent scattering. 

Recentlv |Knez et all J2005h suggested that the ratio in band 
intensity could be due to large amounts of CO2 mixed in with 
the H2O ice. Observ ations reveal that solid CO? is common i n 
many lines of sight JCerakines et al Jll 994 iGih b et al JI2004I) . 
With the Infrared Space Observatory (ISO), the C0 2 stretch- 
ing mode at 4.25 u m (2350 cm" 1 ) was observed toward Tauru s 
background stars JWhittet et alJll998l iNummelin et alJl200ll) . 
More recently the launch of the Spitzer Space Telescope made 
the CO2 bending mode at 15 fim (670 cm -1 ) available for obser- 
vations and t he band has been detected tow ards several back- 
ground stars (Bergin et al. 2005; Knez et al. 2005). The average 
abundance with respect to H2O ice towards the Taurus sources 
is 20%, but up to 37% has been observed (Knez et al., in prep.). 
Toward several pr otostars up to 35% of CO? compared to H?Q 
has been observed ( Nummeli n et al.l200lUBoogert et ali20 04 ), 
making CO2 one of the most abundant ices after H2O. 

In a previous study, using the H2O column density from 
the 3.08 ^m ban d and laboratory spectra of pure H2O ice, 
iKnez etalJJ2005l) determined that the H2O bending mode con- 
tributes 77% and 69% to the observed 6.0 fim absorption fea- 
tures towards Elias 16 and CK 2, respectively. Using a com- 
bination of laboratory spectra of two H20:C02 mixtures, 1 : 1 
and 10:1 respectively, they showed that 85% to 100% of the 
6.0 fim band strength can be explained by H2O. This is due 
to the smaller ratio between the stretching and bending mode 
strengths seen in an unpublished H20:C02 1:1 spectrum by 
Ehrenfreund (private communication). This combination of 
spectra, together with a water-poor mixture spectrum, also fits 
the CO2 profile w ell, with approxim ately 80% of the CO2 in 
the water rich ice jKnez et alJliool . Hence it is not unlikely 
that a significant part of the H2O ice is in H20:C02 ice mix- 
tures close to 1:1 in many lines of sight, even if the average 
abundances of H2O and CO2 is closer to 3: 1 . 

The question that prompted this study is whether there is a 
change in H2O band profiles and relative H2O band strengths in 
ice mixtures compared to pure H2O ice, for H20:C02 ice mix- 
tures that are both astrophysically relevant and contain enough 
H2O to observe with e.g. Spitzer. Such information in turn, 
is a prerequisite for determining whether additional species 
contribute to the 6 fim band. In this work we present a sys- 
tematic study of the infrared properties of H2O absorptions in 
H20:C02 ice mixtures around 1:1, for the temperature range of 



15 to 135 K, in order to constrain the effect that CO2 has on the 
shapes and relative band strengths of the H2O bands. 



2. Previous laboratory data 

Two previous studies have reported changes in the H2O 
bands in H?Q: CO? ice mixtures comp ared to pure H2O ices 
( Ha gen et al]l983l|Schmitt et alJl9 89t. In both cases, only one 
isolated H20:C02 mixture, 1:2 and 10:1 respectively, was in- 
vestigated. In both the 1:2 and 10:1 H20:C02 mixture spectra 
a new H2O band appears around 2.74 fim, the relative strength 
of the bending mode increases and all the band profiles change 
compared to a pure H2O ice spectrum. No attempts were made 
to quantify these changes. 

A number of later laboratory studies focused on H20:C02 
ice mixtures as well, but to our knowledge none of them sys- 
tematically studied the impact of CO2 on water ice and none 
of them has reported on changes in the H2O band profiles and 
band strengths due to CO2. This is not surprising, since most of 
the studies focused on mixtures with H?Q as the dom inant ice 
component ( Hudgi ns et alll993tlBernstein et all2005l) . Only a 
few laboratory spec tra of H?Q:CO? mixtures close to 1:1 ex- 
ist in the literature JGerakines et alJll995t IPalumbo~& Baratta 
200(i lGerakines et alJ2005h The effect of CO? onH ?Q ice fea- 
tures is only considered bv iGerakines et all J2005I) . who con- 
cluded from the spectra of a H20:C02 1 .6: 1 ice mixture that the 
relative H2O band intensities were not affected by high concen- 
trations of CO2. Nevertheless, the 3.08/6.0 fim band strength 
ratio can be calculated from their reported integrated intensi- 
ties and reveals a drop of ~30% compared to pure H2O. 

It is well known from matrix isolation experiments that the 
different H2O bands have dif ferent relative band st rengths de- 
pendent on H2O cluster size ( van Thiel et al. I ll957h . The band 
strength of the bending mode is much less affected by the H2O 
cluster size than the stretching mode, when the H2O cluster 
size is changed by varying the H2O concentration in a N2 
matrix. The former band strength drops by 20% when the 
matrix to H2O ratio is decreased by one order of magnitude 
while the stretching band strength increases by a factor of 10. 
The main intensity contribution in both bands comes from the 
monomer peak when the H2O concentration is low and from 
a cluster mixture when the conce ntration is high (high con- 
centration meaning N2:H20 1:10. Ivan Thiel et aflll957l) . The 
position of the bending modes is approximately the same for 
all cluster sizes, while the major stretching mode peak shifts 
from 3400 cm at high H2O concentrations to the position of 
the free-OH stretch at low concentrations (~3690 cm or ~2.7 
fim). Finally, the relative band strengths are also affected by 
the type of matrix used, e.g. noble gas or nitrogen or oxygen. 
This has been reported also for H? Q in astronomically relevant 
matrices ( Ehrenfreund et al. 1996). It is important to note that 
none of these matrices forms hydrogen bonds with H2O, thus 
no conclusions about the band strengths and band profiles of 
H2O can be drawn from these experiments concerning H2O in 
a hydrogen bonding matrix. Another astrophysically relevant 
molecule, CO, does probabl y form hydrogen b onds with H2O 
in amorphous ice mixtures JSchmitt et alJl9 89). Matrix exper- 
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Table 1. Ice mixtures studied in this work 



Wavelength / /j.m 
4 5 6 



Composition 


H 2 (L) fl 


C0 2 (L) 


Total ice expo 


pure H 2 




3000 





3000 


pure CO2 







3000 


3000 


H 2 0:C0 2 1 


0.25 


3000 


750 


3750 


H 2 0:C0 2 1 


0.5 


3000 


1500 


4500 


H 2 0:C0 2 1 


1 


3000 


3000 


6000 


H 2 0:C0 2 1 


2 


3000 


6000 


9000 


H 2 0:C0 2 1 


4 


3000 


12000 


15000 


H 2 0:C0 2 1 


1 


10000 


10000 


20000 


H 2 0:C0 2 1 


1 


1000 


1000 


2000 



" 1 L (Langmuir) = 1 x 10 6 Torr s ~ 1 monolayer of molecules 

Table 2. The measured peak positions and the integration 
bounds in cirT 1 (jjm) used to compute the integrated intensi- 
ties of the H2O peaks 







Integration bounds 


H 2 bands 


Peak 


Lower 


Upper 


libration 


780(12.8) 


500 (20.0) 


1100(9.09)" 


bend 


1655 (6.04) 


1100(9.09) 


1900 (5.26) 


bulk stretch 


3279 (3.05) 


3000 (3.33) 


3600 (2.78) 


free OH stretch 


3661 (2.73) 


3600 (2.78) 


3730 (2.68) 



" The C0 2 bending mode is excluded by explicitly taking into ac- 
count only half of the libration mode profile (see text) 



iments have also sh own that CO form s weak hydrogen bonds, 
while CO2 does not dTso & Leell985l) . 



3. Experiment and data analysis 

3.1. Experiment 

All experiments were conducted in a high vacuum (HV) c ham- 
ber described in detail elsewhere (Geraki nes et aljfl995l) at a 
base pressures below 1.3xl0~ 6 Torr at room temperature. Ices 
of C 18 C<2 (Praxair 97% purity) and H2O (deionized and further 
purified through sequential freeze-thawing in a vacuum mani- 
fold) were grown on a Csl window, precooled to 15 K (45 K 
for one specific experiment), via effusive dosing at a growth 
rate of ~10 16 molecules cirT 2 s roughly along the surface 
normal. C 18 C»2 was used instead of the main isotopologue of 
CO2 to minimize overlap between H2O and CO2 spectral fea- 
tures. Transmission Fourier transform infrared spectra of the 
ice systems were recorded between 4000 - 400 cirT 1 (2.5 - 25 
fj.m) at a spectral resolution of 2 cirT 1 at fixed temperatures be- 
tween 15 and 135 K, using a total of 256 scans. Background 
spectra were acquired prior to deposition for each experiment, 
at the same resolution and number of scans, and automatically 
subtracted from the recorded ice spectra. 

The pure ices were grown in situ from C 18 C<2 and H2O 
gas bulbs that were filled to a total pressure of 10 mbar, pre- 
pared in a glass-vacuum manifold at a base-pressure of ~1Q 4 
mbar. Mixed ices were made by dosing gas from premixed 
H20:C I8 C>2 bulbs, also prepared in the glass-vacuum manifold. 



0.20 



C 18 stretch 



H„0 free OH stretch 

H bulk stretch 




10 



15 20 



C 1B bend 



H bend 

2 

H libration 




4000 



3000 2000 1000 

Wavenumbers / cm -1 



Fig. 1. Spectra over the 4000-400 cirT 1 range of the H 2 0:C0 2 
1:0.5 ice mixture at 15 K. The assignments of the major peaks 
are indicated. The double headed arrows indicate the integra- 
tion ranges for determining the intensities of the various bands; 
for the libration mode, only half the range is indicated (see 
text). Note that all CO2 frequencies are shifted compared to 
those observed in space, since in the present experiments C 18 C<2 
has been used. 

The ice mixtures studied here are summarized in Tabled 
The relative concentrations for the H20:C 18 C>2 mixtures range 
from 20 to 80% CO2, where the % are relative to the whole 
ice, i.e. 20% CO2 is equivalent to a H20:CC»2 4:1 mixture. The 
H2O ice exposure was always kept at 3000 L (1 L (Langmuir) 
= 1 x 10~ 6 Torr s » 1 monolayer of molecules, assuming 10 15 
molecules cm 2 and a sticking probability of 1), except in two 
1:1 control experiments with 1000 and 10000 L H2O exposure, 
respectively. The nomenclature adopted is as follows, A:B de- 
notes a mixture with A parts of H2O and B parts of C 18 C<2. 

3.2. Data analysis 

The acquired spectra were first reduced to flatten the baseline 
by fitting a second order polynomial to the same 5 points, cho- 
sen by visual inspection well away from any features. This 
was done to facilitate comparison between laboratory and as- 
tronomical spectra. The curved baseline in the raw data is a 
real feature which is due to scattering within the ice. Flattening 
the spectra may hence distort profiles, but in our case no 
such distortion was noted. All reduced spectra can be found 
at http://www.strw.leidenuniv.nl/~lab/databases 

Figure^shows an overview spectrum of a H20:C02 1:0.5 
mixture. The use of the C 18 C<2 isotope changes the position 
of the CO2 bands compared to the main isotopologue, but 
does not affect the conclusions of this study. In addition to the 
C 18 C»2 bending and stretching modes, there are two weak, nar- 
row C 18 C<2 overtone bands in the spectrum, at 3515 and 3625 
cm" 1 respectively, that overlap with the H2O stretching bands. 
To calculate the integrated intensities of the H2O bands, the 
H2O peak intensities were integrated over the same wavenum- 
ber range for all spectra (see Tableland Fig.^. Where narrow 
CO2 peaks overlapped with the H2O bands, a Gaussian was 
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fitted to the CO2 peak and subtracted from the spectrum. A 
Gaussian could not be fitted to the CO2 bending mode, which 
overlaps with the H2O libration mode. The band strength of 
the H2O libration band was instead calculated by doubling the 
band strength of the high frequency half of the band. The peak 
frequency of the H2O libration mode varies in the different mix- 
tures and since the lower frequency integration bound was set 
to the peak frequency, the lower bound is somewhat different 
for different mixtures. For pure H2O ice, this procedure was 
found to be accurate to within 3% compared to integrating the 
whole profile. To further test the induced error of this approach 
the band strengths were also calculated by simply subtracting 
the integrated area of the CO2 bending mode from the total 
integrated area. The difference was less than 5%. The large in- 
terval required for the H2O bending mode is due to its sub- 
stanti al low frequency w ing attributed to the librational over- 
tone JPevlinet all2 001). which in some spectra contains more 
than 50% of the integrated peak intensity. The 2.73 fim band is 
due to free-OH stretches and its assignment to H2O monomers, 
dimers and small multimers is justified in the discussion part of 
this paper. 

The substantial change in the H2O band intensities between 
different mixtures reported here, prevents an independent check 
of the relative amount of H2O deposited onto the surface. From 
three repeated experiments (of the 2:1 mixture) the standard 
deviations of the integrated H2O peak intensities are estimated 
to be less than 10%. These experiments were carried out in the 
beginning and at the end of an experimental series and for sep- 
arately prepared mixtures. The standard deviation hence con- 
tains the error from mixing, absolute ice exposure, measuring 
errors and changes in the experiment over time. Additional er- 
rors may arise from the flattening of the spectra; the difference 
between integrated intensities for the raw and reduced spectra 
is the largest for the 6.0 /mi H2O band at low H2O concentra- 
tions; up to 40% for the 1 :4 mixture. For the astrophysically rel- 
evant mixtures this uncertainty is only 1 - 5%, however. Some 
systematic errors due to the mixing procedure cannot be ex- 
cluded, but are difficult to quantify. Taking all error sources into 
account, we estimate that the relative band strengths are accu- 
rate within ~10% for the astrophysically relevant ice mixtures, 
i.e. ice mixtures with less or equal amounts of CO2 compared 
to H 2 0. 

The band strengths, A, of the three H2O bands present in 
pure H2O ice were estimated for all mixtures using the mea- 
sured band strengths for pure H2O ice at 14 K by Gerakines et 
al. (1995). The band strengths for the pure H2O bands were 
thus set to 2.0 x 10" 16 for the H 2 stretching mode, 1.2 x 
10" 17 for the bending mode and 3.1 x 10" 17 cm molecule" 1 
for the libration band. The band strengths of these bands in the 
H20:C02 mixtures were calculated by scaling each integrated 
intensity by the band strength of the pure H2O ice band over 
the pure H2O band integrated intensity: 

_ ^band 

A H?0:C0 2 =l:x = I / H 2 0:C0 2 = l:xX-^ (1) 

Jhand land ' H 2 

where o co =i x ls tne calculated strength of each H2O band 
in the l:x H20:C02 mixture, Jj^ 7h 2 0:C0 2 =1:x the measured in- 



tegrated intensity of the same band, Aj^ the known strength 

of the pure H2O band and / H ,o the integrated intensity of 
the pure H2O band. The band strengths of the free OH stretch 
in the different mixtures were scaled to the band strength of 
the bulk stretching mode in pure H2O ice. The ratios of our 
measured integrated intensities for pure H2O at 15 K coincided 
with those of Gerakines et al. within 10%. 

The calculated strengths of all H2O bands were plotted as 
a function of CO2 concentration in the ice mixture and fitted 
by linear models as described in detail in section 4.1. These 
models are more accurate in predicting the band strengths for 
H2O in a certain mixture than individual measurements, since 
they are derived from all experiments and hence the random 
errors are averaged out. Our model predictions agreed well 
with previously published isolated measurements. A spectrum 
by Schutte (Leiden Molecular Database) of H 2 0:C0 2 1:1.25 
has a stretching to bending peak ratio which lies within 5% 
of the value predicted by our model fit. In addition, the ratio 
of stretching to bending modes in the H20:C02 1.6:1 mixture 
from Gerakin es et alJd2005l) lies within 8% of our model value, 
further corroborating the results presented in this paper. 

4. Results 

4.1. Changes in H 2 band strengths and profiles with 
mixture composition 

Figure |3 shows that the H2O spectra at 15 K undergo two sig- 
nificant changes as the amount of CO2 is increased from to 
80%; the profiles of the H2O bending and stretching bands 
change dramatically and all H2O integrated peak intensities 
change systematically with varying CO2 concentration. 

Of all modes, the profile the H2O bending mode is most 
affected by the CO2. In pure H2O the bending mode consists 
of a broad peak centered at 1661 cm" 1 (6.02 jjm). The bending 
mode overlaps with the librational overtone, which shows up in 
the spectra as a low-frequency wing. To simplify comparison 
with astronomical observations, the wing is generally treated 
as a part of the bending mode when calculating band strengths. 
The band strength of the pure bending mode is smaller for the 
mixtures that have strong librational modes and hence strong 
librational overtones. In the spectrum of the 1:4 mixture, two 
narrow peaks are observed in the H2O bending region at 1609 
and 1634 cm" 1 , instead of the broad peak and wing in the 
pure H2O spectrum. The profiles of the bending mode in the 
other mixtures appear to be a composition of the peaks in the 
pure H2O spectra and the H20:C02 1:4 mixture. To test this, 
Gaussian profiles were fitted to the bending peak of the pure 
H2O ice, disregarding the wing, and to the two peaks in the 
1 :4 mixture (Lorentzian profiles were also attempted but were 
impossible to match to the ice peaks). The IDL simplex opti- 
mization routine was used to fit two Gaussian peaks simulta- 
neously to the band in the 1:4 spectrum. Figure shows that 
the bending mode in all other mixtures can be separated into 
these three peaks; the slight misfit on the blue wing in all mix- 
tures is corrected for by an additional, small Gaussian centered 
at 1685 cm 1 with FWHM (Full Width at Half Maximum) 35 
cm" 1 . The position of the 1634 cm" 1 peak is shifted with CO2 
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2.70 



2.75 



2.8 



Wavelength / /xm 
3.0 3.2 6 



6.5 



7.5 



12 



14 16 18 



0.06 
0.00 

0.06 
0.00 

0.06 
0.00 

0.06 
0.00 

0.06 
0.00 

0.06 
0.00 



. a) free OH stretch x 4. 
"HO 

2 


1 1 1 

. b) bulk stretch 


1 1 1 — 

c) bend x 5 


1 1 1 1 

d) libration x 5 


'*<~~i*\i i-r-i 1 

- 1:0.25 /v^— \ 




— 1 


1 1 


— 1 1 H" 

. ^~4J 

1 1 1— r 


r — 


' 1 1 

- 1:0.5 




1 


1 1 

1 = 


■===^4- 1 I-- 




1 1 

^^^^^^^ 




— -^1 

— 


^~ 1 


















: -f)k 




: J 


I 




Vi! 



3700 



3650 



3600 3400 3200 1700 1600 1500 

Wavenumbers / cm -1 



900 800 700 600 



Fig. 2. Spectra over the 4000-400 crrT 1 range at 15 K of the four H2O peaks in H20:C02 l:x ice mixtures for different CO2 
concentrations. Each column contains one of the H2O modes and each row one of the H20:C02 mixtures: a) an expanded view 
of the H2O free OH stretch, b) the H2O bulk and free-OH stretch, c) the H2O bend and d) H2O libration. The intensities of the 
libration, bend and free-OH stretch bands have been scaled, with a scaling factor indicated in the first row. The H2O exposure 
was kept constant in all experime nts. The narrow p eaks around 3500 and 3630 crrT 1 have been previously assigned to CO2 
combination modes (Sandford & Allamandola 1990). The sharp feature around 700 cirT 1 is an experimental artifact. 



concentration (2 to 8 cirT 1 for the different mixtures), but the 
FWHMs of the Gaussians were always kept constant. The other 
peak positions were also kept constant. The positions, widths 
and intensities of the Gaussians used to fit the bending mode 
for each mixture are listed in Tabled 

The same bulk stretching peak apparent in the pure H2O 
spectra is visible in all mixture spectra (Fig. [5] column b). 
As more CO2 is mixed in, it acquires more of a double peak 
structure, but the total width of the band remains the same. In 
pure H2O ice the peak at 3696 cirT 1 (2.71 fim) attribut ed to 
the fr ee/dangling OH stretch was not detected ( Rowland et al. 
Il99ll) . In contrast, Fig. |2 column a, shows that a free OH 
H2O stretching band, centered around 3660 cm -1 (2.73 jum), at- 
tributed to the stretches of small H2O clusters, is clearly present 
in all ice mixtures with CO2 and is comparable in intensity to 
the bulk stretching band at high concentrations of CO2. While 
the intensity of the H2O libration mode clearly drops as more 
CO2 is mixed in, it is difficult to tell whether the profile of the 
libration band is affected by the presence of CO2, because of 
the overlap with the CO2 bending mode. 



Of all the H2O bands, only the free OH stretching band 
grows in strength as the amount of CO2 in the ice increases. The 
band strengths of all other H2O ice bands drop with increased 
concentration of CO2, as illustrated in Fig. |4j where all band 
strengths have been scaled to that of the H2O stretch in pure 
H2O ice. The band strength of the bulk stretching band and the 
libration band are strongly dependent on the CO2 concentra- 
tion, while the intensity of the bending mode is less affected, 
whether the librational overtone is included or not, in calculat- 
ing its integrated intensity. The relationship between the inte- 
grated peak intensities and CO2 concentration is well described 
by linear models within this experimental domain and the data 
can be fitted with a typical squared correlation coefficient R 2 = 
0.98 (Table @}. The correlations for the band strengths of the 
bending mode are considerably less, R 2 = 0.81 and 0.89 at 15 
K and 45 K, respectively. This is due to the relatively small 
change in band strength with CO2 concentration. 

In addition, the relative strengths of the various components 
of the bending modes were calculated from the previously fit- 
ted Gaussians. The integrated areas relative to the pure H2O 
bending mode are listed in Tableland plotted in Fig. [5] to- 
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Fig. 3. The three Gaussian components present in the H2O 
bending mode in mixtures with CO2 disregarding the low fre- 
quency wing. In the upper row, first column the pure H2O bend- 
ing mode consists of a wing and a peak, where the peak was 
fitted with a Gaussian at 1661 cirT 1 . In the lower row, second 
column, the H2O bending mode in the 1 :4 mixture can be sep- 
arated into two narrow Gaussians at 1609 and 1634 cm -1 . The 
four other H20:C02 mixtures can be separated into these three 
components derived from the pure and 1 :4 mixed ice. A small 
additional Gaussian centered at 1685 cirT 1 accounts for some 
high frequency excess. The dotted lines indicate the individual 
Gaussian components and the dashed lines their sum. 
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Fig. 4. The band strengths of the four main H2O peaks in 
H20:C02 ice mixtures relative to that of the H2O stretch in 
pure H2O ice at 15 K. The band strengths of the libration, bend 
and free OH stretch modes have been multiplied by a factor of 
10 to facilitate display. The estimated average error bar of the 
relative band strengths is shown in the upper right corner. 

gether with the area of the residual wing. The ratio between 
the 1661 and 1634 cirT 1 peaks is plotted and its significance in 
astrophysical applications is discussed in section 5.2. Table [3] 
also contains the total strength of the bending mode, excluding 
the wing. Due to the exclusion of the wing the drop in band 
strength with CO2 concentration is smaller than in Fig.0] 
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Fig. 5. The integrated area of the three Gaussian components 
present in the H2O bending mode in mixtures with CO2 and 
of the low frequency wing attributed to a librational overtone. 
The ratio between the 1634 and 1661 cirT 1 peaks is plotted as 
well. The estimated error, shown as an error bar in the top right 
corner, is mainly due to the error in the total band strength and 
not to the Gaussian fit. 

Table 5. Ratios between the 1634 and 1661 cirT 1 components 
for all mixtures at 15, 45 and 75 K 



Ratio at given temperature 
Composition 15 K 45 K 75 K 



pure H2O 













H 2 0:C0 2 


1:0.25 


0.34 


0.12 


0.034 


H 2 0:C0 2 


1:0.5 


0.55 


0.28 


0.11 


H 2 0:C0 2 


1:1 


1.1 


0.49 


0.11 


H 2 0:C0 2 


1:2 


2.0 


0.73 


0.17 


H 2 0:C0 2 


1:4 


>10 


1.6 


0.31 



4.2. Temperature dependence 

The profiles of all H2O bands in the H20:C02 mixtures change 
systematically with temperature to become more similar to the 
bands of pure H2O ice as the temperature is increased (Fig. 
HJ. This is especially apparent for the bending mode, where the 
narrow bands lose intensity and the broad band, associated with 
pure H2O ice, gains intensity as the temperature rises. This is 
also seen in Fig.0where the bending mode of the 1 :0.5 mixture 
at 15, 45 and 75 K has been decomposed into the same compo- 
nents derived from the pure H2O and 1:4 mixture at 15 K. As 
the temperature increases the shapes of the components remain 
the same, but the ratio between the thin components and the 
pure H2O component decreases. The ratios for all mixtures at 
15, 45 and 75 K are listed in Table|5] 

The strengths of all H2O bands increase with temperature, 
except for the free OH stretching band which decreases and 
disappears completely above 90 K. Figure [3] shows that the 
bulk stretching band increases monotonically in band strength, 
while the bending and libration bands display jumps and lo- 
cal minima. These jumps are most pronounced around the CO2 
desorption temperature (-90 K). As in the case of concentra- 
tion dependency, the band strengths of the stretching and li- 
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Table 3. The defining parameters - peak position, FWHM (full width at half maximum) and peak height - of the fitted Gaussian 
components of the H2O bending mode. In addition, the integrated area for each peak and their total sum are listed. The area of 
the libration overtone wing is also listed, but not included in the total area. 



Ice mixture 


Peak position 


FWHM 


Relative 


Relative 


Total integrated area 




[crrT 1 (jum)] 


[cm 1 (yum)] 


peak height 


integrated area 


excl. wing 


Pure H 2 


1661 (6.02) 


130 (0.49) 


1 


1 


1 




wing 






0.54 




H 2 0:C0 2 1:0.25 


1609 (6.22) 


15 (0.059) 


0.056 


0.0067 


1.29 




1642 (6.09) 


49 (0.19) 


0.93 


0.31 






1661(6.02) 


130 (0.49) 


0.91 


0.91 






1685 (5.93) 


35 (0.12) 


0.22 


0.055 






wing 






0.60 




H 2 0:C0 2 1:0.5 


1609 (6.22) 


15 (0.059) 


0.14 


0.017 


1.09 




1642 (6.09) 


49 (0.19) 


1.1 


0.36 






1661(6.02) 


130 (0.49) 


0.65 


0.65 






1685 (5.93) 


35 (0.12) 


0.22 


0.055 






wing 






0.32 




H 2 0:C0 2 1:1 


1609 (6.22) 


15 (0.059) 


0.42 


0.050 


0.96 




1638 (6.11) 


49 (0.19) 


1.4 


0.47 






1661(6.02) 


130 (0.49) 


0.40 


0.42 






1685 (5.93) 


35 (0.12) 


0.22 


0.055 






wing 






0.07 




H 2 0:C0 2 1:2 


1609 (6.22) 


15 (0.059) 


1.3 


0.15 


0.99 




1636 (6.11) 


49 (0.19) 


1.6 


0.54 






1661(6.02) 


130 (0.49) 


0.27 


0.27 






1685 (5.93) 


35 (0.12) 


0.14 


0.034 






wing 






0.02 




H 2 0:C0 2 1:4 


1609 (6.22) 


15 (0.059) 


2.8 


0.34 


0.94 




1634 (6.12) 


49 (0.19) 


1.7 


0.57 






1685 (5.93) 


35 (0.12) 


0.14 


0.034 






wing 






0.01 





Table 4. The linear fit coefficients for the H2O band strengths as functions of CO2 concentrations in %, based on six experiments 
with to 80% CO2. The last two rows show the linear fit to the ratio between the bulk stretching and the bending mode. 







Linear Coefficients 




Peak 


Temperature 


constant 


linear coefficient 


R 2 




[K] 


[10~ 16 cm molecule - '] 


[10~ 19 cm molecule -1 ] 




H 2 libration 


15 


0.32+0.02 


-3.2+0.4 


0.99 




45 


0.42±0.03 


-2.7+0.6 


0.92 


H 2 bend 


15 


0.14±0.01 


-0.5+0.2 


0.81 




45 


0.17±0.01 


-0.6+0.1 


0.89 


H 2 bulk stretch 


15 


2.1+0.1 


-22+2 


0.99 




45 


2.8+0.1 


-21+2 


0.98 


H 2 free OH stretch 


15 





1.62+0.07 


0.99 




45 





1.40+0.05 


0.99 






constant 


linear coefficient 


R 1 


H 2 bulk stretch / bend 


15 


16.7+0.2 


-0.160+0.005 


0.99 




45 


17.0+0.3 


-0.101+0.008 


0.99 



bration bands are more affected by temperature than the band 
strength of the bending mode (including or excluding the libra- 
tion overtone wing). Note that the profile of the bending mode 
is the most affected by both changes in temperature and CO2 
concentration, however. 

Different bands are thus differently affected by an increas- 
ing temperature. In addition, the spectra of ice mixtures with 
higher CO2 concentration are more affected by changes in tem- 
perature than those of ices with less CO2; the influence of CO2 



concentration on the intensity ratio between the bulk stretching 
and the bending modes (including the libration overtone wing) 
is plotted in Fig.|9]for different temperatures between 15 and 
135 K. In this plot the CO2 concentration on the horizontal axis 
is the initial one. Above 90 K most of the CO2 has desorbed, 
which changes the mixture composition. In addition, the ratio 
between the stretching mode and the sum of the bending mode 
components (excluding the wing) at 15 K is plotted for com- 
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Fig. 6. Temperature dependence of the H 2 bands in a 1 :0.5 H 2 0:C0 2 mixture i.e. 67% H 2 and 33% C0 2 . Spectra at 4000-400 
crrT 1 of the four main H 2 peaks in a 1:0.5 ice mixture with C0 2 at temperatures between 15 and 135 K. Each column contains 
one of the H 2 vibrational modes and each row one of the temperatures, a) an expanded view of the H 2 free OH stretch, b) the 
H 2 bulk and free-OH stretch, c) the H 2 bend and d) H 2 libration. The spectra of libration, bending and free OH stretch have 
been scaled, with a scaling factor indicated in the top row. 



parison. The parameters of the linear fit to the 15 and 45 K 
measured points are included in Table|4] 

4.3. Dependence on additional parameters: deposition 
temperature and ice thickness 

When the ice mixture (H 2 0:C0 2 1:0.5) was deposited at 45 
K instead of 15 K the profiles and band strengths of the ice 
bands did not change compared to those found when heating 
the 15 K mixture to 45 K. The profiles remain similar also at 
higher temperatures. As the ices are heated, the increase in H 2 
bulk stretching band strength is significantly smaller for the ice 
deposited at 45 K compared to 15 K. In addition, 2-3 times as 
much C0 2 is retained in the H 2 ice at temperatures above the 
C0 2 desorption temperature of 85-90 K (Fig. [TO}. 

Two control experiments with approximately three times 
more and three times less total ice thickness were run to test for 
changes in peak profiles and relative peak band strengths with 
thickness (not shown here). The relative peak band strengths 
did not change significantly for these experiments (i.e. less than 
the previously estimated experimental uncertainty of 10%) . 



Nevertheless, the profile showed some thickness dependence, 
with the narrow peak in the bending mode being more pro- 
nounced for the experiments with less ice coverage. 

5. Discussion 

5.1. Ice structure 

The changes in the H 2 spectral features in H 2 0:C0 2 ice mix- 
tures, compared to the pure H 2 ice, demonstrate that the 
mixed-in C0 2 does affect the H 2 ice structure. The nature of 
the interaction between H 2 and C0 2 is not obvious. One pos- 
sible scenario is that the C0 2 is spread out in the ice more or 
less uniformly and that the change in the H 2 peak intensities 
and profiles is due H 2 bonding with C0 2 . A second option 
is that instead C0 2 forces the H 2 into small clusters, either 
in the gas-phase or upon arrival at the ice surface, since it is 
known from matrix isolated experiments and calculations that 
H 2 mole cules form much s t ronger bonds with e ach other than 
with C0 2 dTso & Ledl 19851 banten et al.ll2005h . The spectral 
changes would then originate from a different type of H 2 0- 
H 2 interaction rather than H 2 bonding with C0 2 . 
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Fig. 7. The components present in the bending mode in the 
H20:CC>2 1:0.5 mixture at different temperatures. The same 
components that were derived in Fig. 0] have been used to fit 
the bending mode at 15, 45 and 75 K. These components are al- 
ways present regardless of temperature, but the ratios between 
the different components changes, especially the ratio between 
the thin 1634 ctrT 1 and the pure H2O 1661 cirT 1 components 
decreases with temperature. 
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Fig. 9. The ratio between the H2O bulk stretching and bending 
mode band strengths (including the libration overtone wing) for 
different temperatures in H20:C02 ice mixtures with increas- 
ing initial amounts of CO2. The CO2 concentration quoted is 
that of the initial ice mixture, while above 90 K most of the 
CO2 has desorbed, changing the composition substantially. In 
addition the ratio between the stretching band and the bending 
mode components (excluding the libration overtone wing) is 
shown for 15 K (multiplied by 0.5 to facilitate viewing). 
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Fig. 8. The band strengths of the four main H2O peaks in an 
initial 1:0.5 ice mixture with CO2 at temperatures between 15 
and 135 K relative to the band strength of the H2O bulk stretch 
in pure H2O ice at 15 K. The band strengths of the libration, 
bending and free-OH stretching mode have been scaled by 10 
to facilitate viewing. The average error bar for the relative band 
strengths is shown in the upper left corner. 
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Fig. 10. The band strengths of the H2O stretching mode de- 
posited at 15 and 45 K in the 1:0.5 ices, respectively, relative 
to the band strength of the H2O stretch in pure H2O ice. The 
CO2 stretching peak is plotted so the amount of CO2 in the ice 
at different temperatures can be compared. The temperature la- 
bels given in the figure indicate the deposition temperature. 



The large decrease in band strength of the H2O bulk stretch- 
ing mode that is observed when CO2 is mixed into the H2O 
ice indicates that CO2 destroys the bulk hydrogen-bond net- 
work, since the band strength of the stretching mode is weaker 
in s mall clusters co mpared to that in larger clusters and bulk 
ice (Ivan Thiel et aljfl957l) . At the same time a new stretching 
peak appears at a higher frequency, i.e. more and more H2O 
molecules are forced into a looser bound environment, which 
increases the intramolecular O-H bond strength, as the CO2 



concentration increases. The frequency of this band, ~2.7 yum, 
and its distinctness from the bulk stretc hing band make it pos - 
sible to assign it to free OH stretches (Rowl and et alJfl99lT) . 
which in general occur both at surfaces of ice and in clusters. 
This agrees with more and more H2O molecules isolated in 
clusters, and hence separated from the H2O ice hydrogen-bond 
network a s CO? is added, but not w ith a uniform arrangement 
with CO2. Ehrenf reund et alJ i 19961) observed a similar band in 
different ice mixtures with 10% or less H2O. They assigned it 



10 



Karin I. Oberg et al.: Effects of CO2 on H 2 band profiles and band strengths 



to overlapping peaks of H2O monomers, dimers and small mul- 
timers consist ent with our analysis. Similar conclusions were 
also drawn bv lvan Thiel et alJ ll 19571) based on their cluster ex- 
periments. 

Simultaneously with the appearance of the free OH stretch- 
ing band, the original broad band in the bending region drops 
in intensity and two narrow peaks appear. The narrowest of 
the two dominates at high CO2 concentration and is at the ap- 
proximate position of the H2O dimer in argon at 1611 cm -1 
(Avers & PullirJ ll976l) . However, the bending mode peaks of 
H2O monomers, dimers and small clusters ar e not far apart and 
shift in position between different matrices llTso & Ledl 19851 
Ivan Thiel et al.lll957i) . thus rather than being due to only H2O 
dimers it is more likely that this peak is produced by a mixture 
of monomer, dimer and small multimer peaks that are over- 
lappin g. The H7O-CO? dimer has its main peak around 1598 
cm" - 1 dTso&Led fl985) and is not visible, further indicating 
that increasing the amount of CO2 lead to a majority of the 
H? Q clustered with themselv es rather than bonded with CO2. 

lEhrenfreund et alJ dl996l) observed a similar bending mode 
to ours in a H20:N2:02 1:5:5 mixture with a narrow peak at 
1606 cirT 1 and a broader feature at 1630 cirT 1 assigned to H2O 
monomers, dimers and small clusters and bulk H2O ice, respec- 
tively. The fact that the same bending mode is acquired when 
H2O is mixed with other small molecules corroborates our in- 
terpretation that the H2O spectrum is dominated by the inter- 
action between the H2O molecules regardless of mixture com- 
position, as long as the mixed-in molecules cannot form hy- 
drogen bonds, or more generally, cannot form bonds with H2O 
of comparable strength to the H2O-H2O bond. Matrix-isolation 
experiments show that e.g. H2O-CO, another astrophysically 
relevant combinati on, form s hydrogen bonds, while H2O-CO2 
does not JTso & Leell985l) . 

As the ice is heated the multimer peaks disappear quickly 
and all other peaks become more like those of pure H2O. As 
the temperature increases the H2O molecules can reorientate 
on their lattice points to a greater degree and form larger and 
larger hydrogen-bonded clusters and finally entire networks. 
This leads to a very small fraction of free OH bonds compared 
to those taking part in hydrogen bonds and hence only bulk vi- 
brations show up in the spectra at high temperatures. Above 90 
K the desorption of CO2 further allows the hydrogen-bonds to 
reform. The temperature affects the H2O molecules directly as 
well as indirectly through its effect on the mobility and final 
desorption of CO2. 

This scenario generally holds independent of whether the 
ice is deposited at 15 or 45 K. At both temperatures the de- 
posited ice should be amorphous and porous, resulting in some 
CO2 molecules becoming trapped in the H2O ice as the ice is 
heated and the pores collapse. The H2O peak profiles are simi- 
lar at all temperatures, but it is clear that more CO2 is trapped 
inside the H2O ice when the ice is deposited at 45 K compared 
to 15 K. Figure ITU1 shows that the higher fraction of CO2 co- 
incides with a lower intensity of the H2O bulk stretching band, 
indicating that the trapped CO2 makes it more difficult for the 
hydrogen bond network to reform. It is not clear at this stage 
why the ice deposited at 45 K is capable, when heated, of trap- 
ping more CO2 within the matrix, but it is most likely a kinetic 



effect, related to the compactness of the ice. The different ice 
structure then affects the relative rates at which pores collapse 
versus CO ? diffusion and desorption during the heating process 
dCollings et all2003l) . 

All experimental results are hence consistent with a model 
in which the H2O molecules are present in the ice either as bulk 
H2O ice or as small H2O clusters. A higher CO2 concentration 
forces more H2O molecules into the cluster state. Since H2O 
bulk ice spectra are different from H2O cluster spectra, adding 
CO2 to H2O ice will significantly change the shapes and band 
strengths of the different H2O peaks compared to pure H2O 
ice. It is expected that other small molecules that cannot form 
H-bonds (or only very weak ones), e.g. O2 and N2, will affect 
H ?Q similarly to CO?, as in dicated by the matrix experiments 
of E hrenfreund et alJdl996l) . Our group is currently conducting 
a systematic series of experiment investigating the effects of 
N2, O2 and CO on H2O spectral features in ice mixtures. 

5.2. Astrophysical implications 

Recent Spitzer observations of the CO2 bending mode at 15 
//m towards background stars indicate that 85 % of th e CO2 col- 
umn density is in a water-rich ice llKnez et all2005t> . To fit all 
H2O and CO2 features consistently, a combination of H20:C02 
10:1 and 1:1 at 10 K (with the ratios of the two mixtures vary- 
ing between 1:0 and 1.3:1 for different sources) was used to fit 
the CO2 profile dKnez et all2005l) . It is hence not unlikely that 
some of the observed H2O ice is in mixtures with close to equal 
concentrations of CO2 and H2O. 

This study shows that CO2 in H20:C02 ice mixtures affects 
the H2O spectrum in three observable ways: 

1 . a free OH stretch peak appears around 2.73 fim, 

2. the profile of the H2O bending mode changes dramatically, 

3. the H2O stretching / bending mode strengths decrease lin- 
early with CO2 concentration. 

When the free OH stretching region at 2.7 yum is observable 
the presence/absence of that band provides a stringent upper 
limit on the amount of CO2 mixed into the H2O ice. This region 
was covered by ISO for high mass YSOs and will be observed 
by Akari (ASTRO-F). It is not covered by Spitzer and cannot 
be observed from the ground. Hence it is not possible to use 
this region of the spectrum to constrain the CO2 concentration 
in H2O ice in most lines of sight. 

Instead we will here use only spectral changes (2) and (3) 
to analyze the Spitzer and Keck spectra of the YSO B5TRS1 
after the silicate absorption has been removed dBoogert et alJ 
2004). We first compare the B5TRS1 spectrum with laboratory 
spectra of pure H2O, CO2 ices at 15 K for reference. Then we 
deduce the maximum amount of H2O that can be present in a 
H20:C02 1 : 1 mixture without a change in peak position for the 
H2O bending mode and then we present the best fit achieved by 
simultaneously comparing the H2O stretching and bending re- 
gions and the CO2 bending mode to laboratory spectra of pure 
H 2 0, CO2 and H 2 0:C02 mixtures at 15 K. An equally good 
fit is achieved using laboratory spectra acquired at 45 K, which 
is also shown. In addition, we show that the average H20:C02 
mixture in the line of sight can be deduced by fitting the 6.0 
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Fig. 11. Spectrum of the YSO B5 IRS 1 (solid line) compared with laboratory spectra (dashed) of pure H2O and CO2 at 15 K (first 
column), a combination of pure H2O and the H20:C02 1:0.5 ice mixture at 15 K (second column) and 45 K (fourth column) 
and the derived bending mode components (third column). The silicate feature has been subtracted from the B5:IRS1 spectrum 
in all plots. In the upper row the spectrum has been fitted with respect to the 3 /urn H2O stretching band. The effect on the H2O 
bending mode region is shown in the second row. In the third row the H2O (mixture) spectra have been subtracted from the star 
spectrum and the residual is compared with a pure HCOOH spectrum at 15 K. The last row shows the fit to the CO2 bending 
mode using the ice mixture derived from the H2O modes. The CO2 peak has been shifted in the laboratory spectra to account for 
the shift in position between C 16 C>2 and C 18 C>2 used in this study. In the fourth column the spectrum of pure HCOOH was fitted 
together with the bending mode components to the 6.0 fim band, to achieve an optimal fit while also being consistent with the 
band strength of the H2O stretching mode. The individual components are plotted with a gray line. 



Table 6. Column densities derived from the different fits to the B5:IRS1 spectra 



Composition 


N(H 2 O tota ,)' 1 


N(H 2 O pure )/ 


N(H 2 Cwd)/ 


N(C0 2 tota) )/ 


N(C0 2 pure )/ 


N(C0 2mixed )/ 




(cm~ 2 ) 


N(H 2 O tota i) 


N(H 2 O tota i) 


N(H 2 O total ) 


N(H 2 O tota i) 


N(H 2 O tota i) 


pure H 2 at 15 K" 


2.1xl0 18 


1 





0.32 


0.32 





H 2 0:(H 2 0:C0 2 1:0.5) 


2.3xl0 18 


0.75 


0.25 


0.23 


0.10 


0.13 


= 3:1 at 15 K c 














H 2 0:C0 2 ~90:10 rf 


2.2xl0 18 





1 






0.08-0.14 


assuming 15 K 














H 2 0:(H 2 0:C0 2 1:0.5) 


1.7X10 18 


0.5 


0.5 


0.20 





0.20 



= 1:1 at 45 K f 



" from H 2 stretching mode 

* first column in Fig. II II 

c second column in Fig. II II preferred fit at 15 K 

d third column in Fig. ll II 

" fourth column in Fig. ll II preferred fit at 45 K 
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yum region with the bending mode components in Fig. [3] while 
keeping the total H2O bending mode integrated intensity con- 
sistent with that of the H2O stretching mode. For all cases we 
also attempt to fit the excess in the H2O bending region with 
HCOOH. The presence of HCOOH is strongly indicated by 
absorption features at 7.25 and 8.2 yum (8.2 yum not shown). 
HCOOH has also been suggested before as th e most li kely 
candi date for the blue wing of the 6.0 fim band JSchutte et alJ 
1996). 

In the analysis of B5:IRS1 we did not take into account 
the effect of grain shapes on the absorption profiles. The effect 
on the H2O bending profile should be small, however, due to 
that the imaginary part of the optical constant (k) for the H2O 
bending profile is small. When k is small it induces (through 
Kramers-Rronig relation) only small fluctuations of n around 
the central peak wave number so that particle scattering will not 
induce a significant change in band shape compared to trans- 
mission spectra. 

In Fig. the laboratory spectrum of pure H2O is scaled 
to the 3.0 yum stretching band (upper row, first column). The 
fit in the bending mode region is shown in the second row and 
the spectrum after subtraction of the pure H2O spectrum in the 
third row. A pure HCOOH spectrum was then scaled to the 
residual. The lower row shows the fit of the CO2 bending mode 
when only a pure CO2 spectrum is used. It is found that the fit 
between the pure laboratory ice spectrum and observed spec- 
trum is poor. We subsequently tested the maximum amount of 
H2O that can be mixed in the 1:1 mixture with CO2 without 
visibly changing the H2O bending mode profile. If 6% of the 
total H2O column density is in a 1:1 mixture there is a clear 
change in peak position of the bending mode so 5% was set as 
an upper limit. 

In the second column a combination of pure H2O and the 
H 2 0:C0 2 1:0.5 mixture at 15 K was used to fit the H 2 stretch- 
ing and bending regions simultaneously. A good fit is also pos- 
sible using the H20:C02 1:1 and 1:0.25 mixture. The best fit 
is achieved using 3 parts pure H2O ice and 1 part 1:0.5 mix- 
ture. The pure/mixture combination results in a smaller residual 
in the H2O bending region when subtracted from the B5:IRS1 
spectrum, compared to when only a pure H2O spectrum is used. 
In addition, the fit with the spectrum of HCOOH to this residual 
is better. The same 1 :0.5 ice mixture, together with pure CO2, is 
also used to fit the CO2 15.2 yum bending mode. Approximately 
50% of the CO2 is found to be in a H2O rich ice for the best fit. 
A combination of pure and 1 :0.25 mixture also results in a good 
fit for the CO2 bending mode, but a combination of pure ice 
and the 1 : 1 mixture. Hence fitting the H2O and CO2 ice modes 
simultaneously effectively constrains both the abundances and 
different environments of both H2O and CO2 ice. 

The third column shows the best fit to the 6.0 yum region 
using the H2O bending mode components derived previously 
and pure HCOOH. This method of estimating the composition 
of water rich ice has the advantage that no mix and matching 
with different laboratory spectra is needed and the uncertainty 
is easier to estimate. The three components were varied inde- 
pendently, but the wing residual was scaled to the pure H2O 
component. The ratio between the broad pure H2O peak at 6.02 
fj.m (1661 cirT 1 ) and the 6.09 fj.m (1642 cirT 1 ) peak was used 



to constrain the amount of CO2 present by comparing this ra- 
tio to Fig. [5] The total H2O bending mode integrated intensity 
was constrained by the ratio between the stretching and bend- 
ing mode band strengths including the wing (Table |4} and the 
ratio between the broad pure H2O component and the narrow 
H2O cluster bending component (Tabled. The latter compari- 
son constrains the amount of H2O mixed with CO2, which must 
be fed back into the ratio between the stretching and the total 
bending mode strength to calculate a correct column density. A 
good fit was achieved for average H20:C02 mixtures 92:8 to 
86:14, which is consistent with the derived composition from 
the best fit in the second column. When higher resolution spec- 
tra of this region will be available, it should be possible to ana- 
lyze all the different components of the H2O bending mode as 
demonstrated in the laboratory spectrum (Fig. 01. This will fur- 
ther constrain the mixture composition from the ratio between 
the two narrow bending peaks apparent in all mixture spectra. 
In the current spectrum the ratio between the two narrow peaks 
only puts an upper limit of 1 : 1 on the dominant H20:C02 mix- 
ture in the line of sight. The column densities derived from 
the fit of laboratory spectra and bending mode components are 
shown in Table [6] 

All the fits above were done with spectra at 15 K. If the ob- 
served ice is at a higher temperature, the same amount of CO2 
mixed into the ice will result in a H2O bending mode profile 
that is less different from pure H2O ice (Fig. [6ji. The best fit 
at 45 K has an average mixture of H20:C02 ~ 1:0.23 (with 
50% of the H 2 ice in a 1:0.5 mixture and 50% pure H 2 0) and 
then no additional pure CO2 is necessary to fit the CO2 bending 
mode. A similar result is achieved when fitting the H2O bend- 
ing mode components and calculating the composition using 
the values of the component ratios for 45 K instead of 15 K in 
Table[5] More generally the ratio of the components constrains 
a region in the ice mixture-temperature space (initial mixture if 
the temperature is above 90 K) rather than determining an exact 
mixture composition. The fit to B5TRS1 resulted in a ratio of 
the 1634 to 1661 cirT 1 components of 0. 15. According to Table 
[5]this means that the temperaturexomposition is bound by 15 
K:ll+3%, 45 K:23±3% and 75 K:39+5% C0 2 relative to the 
H2O ice abundance. Using the components to fit the bending 
region rather than spectra is advantageous if the temperature is 
unknown, since the ratio between the components sets bounds 
on the temperature and composition simultaneously and fully 
reveals the uncertainties involved. In most sources it is possi- 
ble to independently constrain the temperature, however. For 
B5TRS1 Boogert et al. (2004) have concluded that the bulk of 
the ice must be colder than 50K. For these temperatures the 
ratio of the bending mode components results in an average 
H 2 0:C0 2 mixture of 1 1 to 23% C0 2 with respect to the H 2 
ice. This can be further refined by more detailed radiative trans- 
fer modeling of the dust temperature in the source by fitting the 
spectral energy distribution. In contrast the total abundances of 
H2O and CO2 are only weakly dependent on temperature and 
can be determined within 15% uncertainty. 

Note that the strengths of the stretching and libration bands 
are affected by CO2, which needs to be taken into account when 
deducing the H2O column density from the optical depths of 
the H2O bands if a considerable amount of the H2O is in mix- 
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tures with CO2. In the example shown here the drop in band 
strength of the stretching band is only ~5% compared to the 
pure H2O ice (Table |6jl. In the case where ~50% of the H2O 
ice is in a 1:1 mixtu re with CO?, as ha s been deduced for some 
background stars bv lKnez et alJ ( 120 05), the band strength of the 
H2O stretching mode would drop by as much as ~25% com- 
pared to pure H2O ice and the inferred H2O column density 
would be ~33% higher. 

Whether the results presented here can explain the ratio be- 
tween the H2O stretching and bending modes in all astrophys- 
ical sources has to be determined for each source individually 
because of the drastic change in the bending mode profile that 
accompanies the drop in ratio as the concentration of CO2 is in- 
creased. In the case of B5:IRS1 the H2O and CO2 band shapes 
and intensities are well described by 50% of the CO2 and 25% 
of the H2O ice in a h^OiCCh 1:0.5 ice mixture at 15 K and by 
100% of the C0 2 and 50% of the H 2 ice in a H 2 0:C0 2 1:0.5 
ice mixture at 45 K. These ratios in turn can provide interesting 
constraints on the H2O and CO2 formation routes. 



6. Conclusions 

1. The H2O band profiles and strengths in astrophysically 
plausible mixtures with CO2 are significantly different 
compared to pure H2O ice. The changes in H2O band pro- 
files with CO2 are greatest for the H2O bending mode, 
while the stretching band displays more subtle changes. 

2. The band strengths of all major H2O bands depend lin- 
early on CO2 in the investigated domain of 20% to 80% 
of CO2 in the ice mixture. The band strength ratio between 
the H2O bulk stretching and bending modes is also linearly 
dependent on CO2 concentration in astrophysically relevant 
H20:CC>2 mixtures such that a H20:C02 1:1 mixture has a 
ratio of half that of pure H2O. 

3. The H2O bending mode in H20:CC<2 mixtures can be sep- 
arated into three components due to pure H2O ice, large 
H2O and small H2O clusters, respectively. In addition, a Vi- 
brational overtone overlaps with the bending mode. 

4. The H2O bending profile is visibly affected by as little as 
6% H2O in a H20:C02 1:1 mixture. The composition of 
the H2O rich ice is hence crucial to identify other species 
contributing to the 6 fim band. 

5. A free OH stretching band appears with CO2 concentration 
which can be used to put strict upper limits on CO2 mixed 
with H2O whenever the 2.7 ^m region is observed. 

6. If the ice temperature can be estimated independently, the 
total amount of CO2 and H2O as well as the amount of CO2 
mixed with the H2O ice can be constrained by fitting the 
appropriate laboratory mixture spectra consistently to the 
H 2 and C0 2 bands. 

7. The average H20:C02 ice mixture in any line of sight can 
be efficiently deduced by fitting the H2O bending mode 
with components corresponding to pure and CO2 rich ice. 
Especially if the ice temperature is uncertain, fitting the 
H2O bending mode components to the bending mode re- 
gion together with HCOOH, will effectively constrain the 
composition and temperature simultaneously. 



8. The deduction of mixture composition is substantially de- 
pendent on the assumed dust temperature. The determina- 
tion of total column densities of H2O and CO2 are only 
weakly dependent on the assumed dust temperature, how- 
ever. 
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